Interactions of distinct quadrupolar nematic colloids 
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The effective interaction between spherical colloids in nematic liquid crystals is investigated in the 
framework of the Landau-de Gennes theory. The colloids differ through their interaction with the 
nematic. While both particles induce quadrupolar far-field distortions in the nematic matrix, with 
unlike quadrupole moments, one favours homeotropic and the other degenerate planar anchoring of 
the nematic director. In the strong anchoring regime the colloids with homeotropic anchoring are 
accompanied by an equatorial disclination line defect, known as "Saturn-ring", while the colloids 
with degenerate planar anchoring nucleate a pair of antipodal surface defects, called "Boojums" . In 
the linear (large-distance) regime the colloidal interactions are of the quadrupolar type, where the 
quadrupoles have opposite signs. These are attractive when the colloids are aligned either parallel 
or perpendicular to the far-field director. At short distance, non-linear effects including "direct" 
interactions between defects give rise to a repulsion between the particles, which prevents them from 
touching. This finding supports the stability of nematic colloidal square crystallites the assembly of 
which has been reported recently. 

PACS numbers: 61.30.-v, 61.30.Jf, 82.70.Dd 
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1. INTRODUCTION 

The self-assembly of colloidal particles into structures 
with controlled spatial ordering is of great importance in 
colloid science, with particular interest in the assembly 
of photonic crystals [l[ - artificially produced periodic 
dielectric structures designed to control and manipulate 
light. In this context, a variety of colloidal structures 
assembled in liquid crystal (LC) matrices @, H[ , combined 
with the unique mechanical and electro-optical properties 
of the LC host [J] , have proven to be good candidates for 
the development of colloid crystals with tunable photonic 
properties. 

In conventional colloids, in isotropic fluids, the col- 
loidal particles interact via van der Waals, electrostatic, 
or steric forces. These forces are isotropic, and their 
range does not exceed a few tens of nanometers. By 
contrast, when dispersed in a LC, due to its long-range 
orientational molecular ordering, colloidal particles inter- 
act predominantly through long-range anisotropic forces 
0, [a] ■ The origin of these effective forces is the elastic 
distortions of the LC matrix due to the presence of the 
colloidal particles. The range of the elastic forces is of the 
order of several colloidal diameters. Elastic forces drive 
the particles to self-assemble into linear chains |7HlO| , pe- 
riodic lattices Iipl-fl2l . anisotropic clusters [13[, and cel- 
lular structures [14| . 

A distinctive feature of LC colloids is the presence 



nunos@cii.fc.ul.pt 



1 miko@is.mpg.dc 



of topological defects [l5[, that not only determine the 
symmetry of the long-range colloidal interaction [l6|, [TtJ , 
but also stabilize the ordered aggregates at short range 
, where the elastic interactions are dominated by non- 
linear effects and render the self-assembly of LC colloids 
a challenging theoretical problem. Topological defects in 
nematic liquid crystals (NLC) are nucleated due to the 
mismatch of the global and local (at the colloidal sur- 
faces) boundary conditions leading to frustration of the 
uniform nematic order. 

Small spherical particles (< lyum) imposing 
homeotropic surface anchoring on the nematic di- 
rector stabilize equatorial Saturn-ring defects [l8j . 
which for larger particles may be stabilized through 
confinement or by external electric fields (20| . The 
far-field distortions have quadrupolar symmetry and 
the resulting pairwise colloidal interaction is of the 
quadrupolar type decaying with the distance d between 
the particles as d~ 5 [lq ]. 

Larger particles (> 0(l/im)) with homeotropic anchor- 
ing induce point-like hedgehog defects [2l|, which lead 
to a far-field director of dipolar symmetry and a large 
distance colloidal interaction varying as d~ 3 For 
particles with planar degenerate anchoring, two antipo- 
dal surface defects, known as boojums [21[ are nucleated 
which lead to a far-field director of quadrupolar symme- 
try. Recent numerical calculations reveal that the cores of 
nematic boojums can take three different configurations 
[22l |: single core, double core, or split core. The single- 
core boojum is a point-like index 1 defect with azimuthal 
symmetry, the split-core boojum has two index 1/2 sur- 
face point-like defects connected by a 1/2 bulk disclina- 
tion line. The double-core boojum is an intermediated 
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FIG. 1: Schematic representation of interacting Saturn-ring 
and boojum quadrupolar particles of radius R. d is the inter- 
particle distance and 9 is the polar angle relative to the far- 
field director. 



structure with broken azimuthal symmetry. The far-field 
distortions and the resulting asymptotic pair interaction 
between particles with degenerate planar anchoring are 
of the quadrupolar type |23(, as for Saturn-ring particles. 

Recently, a light sensitive coating of colloidal parti- 
cles was used in order to switch the surface anchoring 
from homeotropic to planar, reversibly J24[, providing 
the means for tuning the colloidal interactions and thus 
controlling the assembly of the colloidal structures. 

Direct assembly of two-dimensional (2D) crystallites of 
quadrupolar dipolar [25| , or dipolar and quadrupo- 
lar [26|, spherical particles has been achieved by using 
laser tweezer techniques. The crystallites are stabilized 
by the presence of topological defects, which provide 
local free energy minima of the order of ~ lOOOfc^T, 
where ks is the Boltzmann constant and T the abso- 
lute temperature. By contrast, interacting quadrupolar 
boojum-particles in three-dimensional systems do not ex- 
hibit short-range repulsive behaviour, and the equilib- 
rium configuration corresponds to close contact or coa- 
lescence of the particles [22|, [23, HI] . 

Recently, the assembly of 2D colloidal crystallites 
of spherical particles with two types of anchoring, 
homeotropic and degenerate planar, dispersed in 5CB 
has been reported [29(. The particles with radius R = 
2.16^im were confined to a cell of thickness 6^m in or- 
der to stabilize the Saturn-ring configuration around the 
particles with homeotropic anchoring. Both types of par- 
ticles generate quadrupolar nematic distortions, but the 
corresponding quadrupolar moments have different signs, 
meaning that the particles attract each other when they 
are aligned parallel or perpendicular to the far-field di- 
rector. This allowed the assembly of 2D square colloidal 
crystallites. In the following we shall adopt the notations 
of Ref. [29I ]: a Saturn-ring quadrupolar particle will be 
denoted by "S", and a boojum quadrupolar particle by 
"P" (planar anchoring). In this article, we present the 
results of a numerical study of the NLC-mediated inter- 
action between S and P particles (see Fig. [T]) for a wide 
range of distances. We focus on the short-distance regime 
that determines the final equilibrium configuration and 



its stability. 

The paper is organized as follows: in the next section 
we introduce the Landau-de Gennes free energy func- 
tional, which we minimize numerically. In Sec. [3] we 
present our results. We show that the S-P interaction 
exhibits two local minima for configurations in which 
particles align either parallel or perpendicular to the far- 
field director, thus allowing the assembly of square lat- 
tices. We show that the two equilibrium configurations 
are separated by a free energy barrier of the order of 
250fcsr(i?/l/^m), where R is the colloidal radius. In 
Sec. 2] we present our conclusions. 

2. LANDAU-DE GENNES THEORY 

Phcnomcnologically, nematic liquid crystals are char- 
acterised by a traceless, symmetric tensor order- 
parameter Q which can, in general, be written as [30l | 

Qij = 2 ( 3n i n i ~ <%) + ~K (hlj ~ miirij) , (1) 

where Q is the uniaxial order parameter, which measures 
the degree of orientational order along the nematic di- 
rector n, and B is the biaxial order parameter, which 
measures the degree of orientational order along the di- 
rections perpendicular to n, characterised by 1 and m. 
The corresponding Landau-de Gennes (LdG) free energy 
functional is [3l[ 

F= [ d 3 x(f b + f e ) + f daf a , (2) 
Jn Jon 

where f b and f e are the bulk and elastic free energy den- 
sities, given by 

ft = a(T) TrQ 2 - 6 TrQ 3 + c (TrQ 2 )) 2 , (3) 

fe = ^(d k Qijf. (4) 

The first integral in Eq.® is over the volume occu- 
pied by the nematic, O, and the second is over the sur- 
faces of the colloidal particles, DO.. The bulk parameter 
a(T) = ao (T — T*) depends linearly on the temperature 
T, with ao a material dependent constant and T* the 
supercooling temperature of the isotropic phase; b and c 
are positive material dependent constants. For a given 
temperature T and in the absence of elastic distortions 
the free energy density in Eq. is minimized for the 
(bulk) degree of orientational order given by 

where r = 24ao(T — T*)c/b 2 is the reduced tempera- 
ture, r controls the stability of the nematic and isotropic 
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phases. For r > 9/8 the nematic phase is unstable, while 
the isotropic phase is unstable for r < 0. The nematic 
and the isotropic phases coexist at r = 1 . For simplicity, 
we consider the one-elastic constant approximation where 
L is related to the Frank-Oseen (FO) elastic constant K 
by K = 9Q 2 L/2 

We consider the interaction between colloidal particles 
with homeotropic and degenerate planar anchoring. We 
assume rigid homeotropic boundary conditions, while the 
planar degenerate anchoring is described by the covariant 
surface anchoring free energy proposed by Fournier and 
Galatola (FG): 



Wi (Qa - Qt) 2 + W 2 (q% - (3Q b /2 



(6) 



Qij = Qij + Qb^-,Qij = (Sik - v t v k )Q kl (Sij - viUj), 
and v is the surface normal. The quadratic term favours 
tangential orientation of the director n, and the quartic 
term guarantees the existence of a minimum for the scalar 
order parameter at the surface equal to its bulk value 
Qb. For simplicity we will assume W\ = W2 = W. The 
equilibrium surface orientation is determined by the free 
energy functional minimum that satisfies the far-field 
boundary conditions. 

Typical values of bulk parameters for 5CB arc ciq = 
0.044 x 10 6 J/Km 3 , b = 0.816 x 10 6 J/m 3 , c = 0.45 x 
10 6 J/m 3 , L = 6x 10" 12 J/m, and T* = 307K. At the 
nematic-isotropic coexistence the bulk correlation length 

1 /2 

£ = (24cL/6 2 ) ~ lOnm determines the order of mag- 
nitude of the spatial extension of the cores of the topo- 
logical defects [32| . Calculations were performed for par- 
ticles of radius R = 0.1 /im~ 10£, in a cubic box of length 
I = 30R; the reduced temperature r ~ 0.16, and the 
anchoring strength W satisfying WQ\RjK ~ 37, which 
corresponds to the strong anchoring regime. 

The Landau-de Gennes free energy Eq.Q is minimized 
by using finite elements methods, with adaptive meshing. 
The surfaces of the spherical particles <9f2 are triangu- 
lated using the open source GNU Triangulated Surface 
Library |33j. Then, the triangulation of the nematic do- 
main Q, is carried out using the Quality Tetrahedral Mesh 
Generator I34| . which supports the adaptive mesh re- 



finement. Linear triangular and tetrahedral elements are 
used in 2D and 3-D, respectively. Generalized Gaussian 
quadrature rules for multiple integrals [35j are used in 
order to evaluate integrals over the elements. In particu- 
lar, for tetrahedra a fully symmetric cubature rule with 
11 points [36[ is used, and integrations over triangles are 
done by using a fully symmetric quadrature rule with 7 
points [37| . The discretized Landau-de Gennes functional 
in then minimized using the INRIA's M1QN3 [38| opti- 
mization routine, which implements a limited memory 
quasi- Newton technique of Nocedal (3|| . 



3. RESULTS 

Lets us assume that the far-field director is parallel 
to the z— axis, and that S, or P particle is placed at the 
origin of the reference frame. For spherical particles, the 
distortions of the director field n(r) are uniaxial, with 
the symmetry group Coo V . At large distances r from the 
particle, the director exhibits small distortions from its 
uniform far-field alignment n(r) ~ {n\, ri2, 1 — 0{n\, 
In the one-elastic constant approximation, the trans- 
verse components (i = 1, 2) satisfy Laplace's equation 
An, = [2j|, and the asymptotic solution for n, can be 
expanded in terms of multipolcs. Due to the symmetry 
requirements, the quadrupolar term is the lowest-order 
term in the expansions [23[, i.e., 
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and the symmetry of the director requires the quadrupole 
moment tensor Qi a p = QSiabsp. In Ref. [23[ an explicit 
expression for Qi a p in terms of the surface integrals in- 
volving ni has been given. Applying then a superpo- 
sition approximation, the pair interacti on p otential be- 
tween quadrupolar particles is obtained j23l l40l| 
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Q1Q2 (9 -90 cos 2 (9 + 105 cos 4 (9), (8) 



where the subscript indicates that the interaction is be- 
tween two uniaxial quadrupoles with quadrupole mo- 
ments Q\iQi- The interaction potential in Eq. @ is 
repulsive at 9 = and 9 = ir/2, and attractive at interme- 
diate orientations. The effective interactions between two 
S, or two P particles have been analysed theoretically be- 
yond the superposition approximation in Refs. (5. l22i. l28l] . 
In both cases, significant deviations from the asymp- 
totic quadrupolar behaviour have been found at short 
distances. 

The interaction between S and P particles was mea- 
sured recently (29|. Asymptotic analysis suggests that S 
particles induce elastic quadrupoles with quadrupole mo- 
ment Qs > 0. By contrast, P particles are characterised 
by quadrupole moments Qp < [23[. The sign of a 
quadrupole moment cannot be obtained from the interac- 
tion of like quadrupoles, S-S or P-P, and only its absolute 
strength can be extracted from such experiments. Differ- 
ent experiments yield a range of values, as presented in 
Table |IJ The discrepancies are most likely due to the dif- 
ferent conditions (e. g. geometrical confinement) under 
which the experiments were performed (4l| . or to the ef- 
fects of thermal fluctuations that are important at large 
distances 27). These factors may limit (severely) the 



range of validity of the asymptotic expression, Eq.©. 

Recent numerical calculations revealed that under cer- 
tain conditions the boojums of P particles have a split 
core structure with broken axial symmetry [22J. In this 
case, the P particle induces a biaxial quadrupole [23[ 



TABLE I: The absolute values of the elastic quadrupole mo- 
ment \Qp\ as reported in different experiments for the case of 
spherical colloidal particles with planar degenerate anchoring. 
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thus adding a correction to the asymptotic result, Eq.©. 
However, since the core size is of the order of the nematic 
correlation length, the effects of the biaxiality should be 
small far away from the particle. Indeed, as we shall 
demonstrate below, our calculations are in good agree- 
ment with Eq. ([8]) indicating that the biaxial contribution 
is small. 

We calculate the two-body effective interaction poten- 
tial (d, 0) [43| which is defined by decomposing the 
excess (over the free energy of the uniform nematic) free 
energy for two particles F2 = F — fb(Qb)Q as follows 

F 2 (d,9) = F^ +F p 1] +F^{d,9). (9) 

F^p is the excess free energy calculated independently 
for an isolated S, P particle, respectively. By definition, 
F^ tends to zero when d — > 00. 

Figure [2] illustrates the S-P interaction potential 
F( 2 \d) for several values of 0. By contrast to the S-S 
and P-P potentials, the S-P interaction is purely attrac- 
tive at 9 = 0, 7r/2, and repulsive at intermediate orien- 
tations, which is consistent with QsQp < 0. As shown 
in Fig. [2 the asymptotic result Eq.® is obeyed accu- 
rately for d > 3R. For smaller distances the superposi- 
tion approximation is qualitatively wrong, indicating the 
significance of non-linear effects. F^ (d, 9) exhibits a lo- 
cal minimum at d± ~ 2.25R, corresponding to 9 = tt/2, 
and a global minimum at d|| ~ 2.16R for relative orienta- 
tion 9 = (see also Fig. 2]). These values are in excellent 
agreement with the experimental results d e * p ~ 2.28R 
and dT p ~ 2.16i? [||. At = tt/4, F (2 1(d) reveals a 
crossover from repulsive, at large d > 2.9R, to attractive 
behaviour, at smaller d < 2.9R. The crossover is driven 
by the structural re-arrangement of the defects (see the 
configurations in Fig. and b) below). This is simi- 
lar to the crossover observed in the S-S interaction [|| at 
= 7r/2, where the Saturn-rings deform from their equa- 
torial position, or for the P-P interaction at = 0, where 
the inner boojums slide on the colloidal surfaces in oppo- 
site directions [22|, [28[ . In the S-P interaction potential, 
however, we observe both the sliding of the inner boo- 
jum and the deformation of the Saturn-ring, as shown in 
Fig. g] below. The configuration at = 7t/4, d ~ 2.14i? 
does not correspond to a local minimum of the S-P in- 
teraction potential. In fact, F^(d,9) in the vicinity of 
that point has a cusp-like shape in the 0— direction (see 
Fig. At very small distances, however, the conflicting 
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FIG. 2: Pair interaction potential as a function of inter- 
particle separation d for different orientations 6 = 0, n/4, n/2 
relative to the far-field director (see FigfTJ . Colored lines rep- 
resent the asymptotic result, Eq.([8j), with QsQp = — 0.18i? 6 . 
The inset shows in more detail the global (9 = 0) and lo- 
cal (0 = 7r/2) minima located at d = du — 2.16R and 
d = d± ~ 2.25-R, respectively. The nematic configurations 
are illustrated in a) for the global minimum at 6 = 0, d = d|| 
and in b) for the local minimum at 9 = n/2,d = d±. The 
iso-surfaces correspond to Q = 0.25, Qb — 0.44, and the short 
black lines represent the director field in the plane y = 0. 



anchoring conditions on the particle surfaces lead to a 
repulsive F^ 2 \d) at any orientation 0. 

This can be seen also from the behaviour of the elas- 
tic force F(d, 0) between the colloidal particles. Fig- 
ure [3] shows the radial F r = —dF^/dd and the polar 
Fg = —(l/d)dF^/d9 components of the force as func- 
tions of d for several values of 0. It is clear that for the 
orientations = and 9 = n/2 the interaction is purely 
attractive at large separations d, and that in all direc- 
tions a strong repulsion appears at very small d, stabi- 
lizing the dimer at a fixed separation d > 2R. Moreover, 
the behaviour of Fg clearly shows that once the particles 
are oriented along = Oor0 = 7r/2 they will approach 
each other in this directions under the action of F r , and 
will stabilize at d = d^ or d = d±, respectively. The ra- 
dial component of the force at = 7r/4 is repulsive for 
large inter-particle distances, and turns into attractive 
as soon as the defects start to rearrange. Although, the 
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FIG. 3: Radial F r and polar Fg components of the force acting 
on the colloidal particles as a function of the inter-particle 
distance d, for several orientations 9 — 0, tt/4, tt/2 relative to 
the far-field director rioo. 



asymptotic result Eq.© predicts a non-zero Fg < at 
9 = 7r/4, i.e., a polar component of the force which tends 
to align the particles parallel to n^, it is only when the 
non-linear elastic effects become important (d < 3i?) that 
Fg increases significantly, and becomes comparable to F r . 
This is in sharp contrast to the P-P interaction force [22j , 
where Fg is always two orders of magnitude smaller than 
F r . 

Figure 0] illustrates the interaction potential (9) 
for several values of the inter-particle distance d. At 
d = 3.0R the interaction potential is a smooth function of 
9. It is clear that for this distance the asymptotic result, 
Eq.©, represented in Fig. fusing the experimental value 
QsQp = — 0.18i? 6 , is qualitatively but not quantitatively 
accurate. As mentioned previously, the asymptotic result 
describes the full LdG results quantitatively at distances 
d > 3R. As the distance between the particles decreases, 
(9) reveals a cusp at some oblique angle 9. This cor- 
responds to a structural transition between two different 
nematic configurations, which "coexist" at the orienta- 
tion corresponding to the cusp. Two representative ne- 
matic configurations are shown at d = 2AR in Fig. 0k) at 
9 ~ 0.55(tt/2) and b) at 9 ~ 0.62(tt/2). When the S and 
P colloidal particles are at the equilibrium configuration 
at 9 = 0, increasing 9 forces the inner boojum to slide in 
the same direction as a result of the repulsion from the 
disclination ring, which deforms upwards from its equa- 
torial position, as depicted in Fig. 0^). The S-P system 
follows the left branch of F^ {9) as 9 increases from zero. 
This configuration becomes metastable at the cusp and 
at a larger value of 9*(d) the repulsion between the in- 
ner boojum and the disclination ring renders it unstable. 
Then, the system jumps to the right branch, and assumes 
a configuration similar to that depicted in Fig. 0b) . On 
the other hand, when the S-P system starts from the lo- 
cal equilibrium configuration at 9 = 7r/2, it follows the 
right branch of F^ as 9 decreases. The disclination ring 
now deforms downwards and pushes one of the boojums 
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FIG. 4: Pair interaction potential f' 2 ' as a function of 
the relative orientation 9 for different inter-particle distances 
d. The black line shows the asymptotic result Eq.([8j) with 
QsQp = — 0.18i? 6 at d = 3R. The snapshots depict the 
nematic configurations obtained from the numerical calcula- 
tions at d = 2AR a) 9 ~ 0.55(tt/2) and b) 9 ~ 0.62(tt/2) 
close to the free-energy cusp. The isosurfaces correspond to 
Q = 0.25, Qb — 0.44, and the short black lines represent the 
director field in the plane y = 0. 



away from the z— axis as depicted in Fig. HJd). 

It is clear from Figs. [5] and 0J that at short-distances 
d < 3R the inter-particle interaction deviates signifi- 
cantly from the asymptotic result described by Eq.©. 
In Fig. [5] we plot the absolute value of the S-P interac- 
tion potential as a function of d for 9 = 0, tt/2. For large 
separations, d > 3R, F^ 2 \d) can be described by a power 
law 



F (2) = f(0)(R/d) aW , 



(10) 



where both a(9) and f{9) are fitting parameters. Figure 
[^indicates that the asymptotic result is valid for d > 3R. 
Thus, in Fig. [5] we fit Eq. (fTU|) to data corresponding to 
d > 3.5R only. For both values of 9 we find a = 4.8 ±0.2, 
which is in good agreement with the asymptotic power 
law oc 1/g? 5 . From these fits we also obtain QsQp/R 6 = 
-0.15 ± 0.03 and QsQp/R 6 = -0.20 ± 0.05 for 9 = 
and 9 = tt/2, respectively. This is in excellent agreement 
with the experimental value QsQp/R 6 = —0.18 reported 
in Ref. |H. 
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such distances the topological defects start to interact 
and re-arrange their positions (see Fig.[3]a), b)), resulting 
in profound changes of the interaction potential (d) 
(see Figured squares). 

Under some conditions, the P particles nucleate boo- 
jums of complex biaxial structure p3 |. This leads to a 
biaxial quadrupolar distortion, since it breaks the axial 
symmetry, and could change the asymptotic uniaxial ex- 
pression given by Eq. ((5]) [2j]. However, since the size of 
the biaxial regions is of the order of the nematic corre- 
lation length, the biaxial perturbation at large distances 
are negligible, which explains the good agreement be- 
tween our numerical results and the asymptotic expres- 
sion for uniaxial quadrupoles. 



FIG. 5: Pair interaction potential as a function of the 
inter-particle distance d for 8 = 0, tt/2. Black lines corre- 
spond to a power law fit F^ = f(9) (R/d) a ^ to the numer- 
ical results. In both cases a = 4.8 ± 0.2. 



4. CONCLUSIONS 

Recently, a directed assembly of spherical particles 
dispersed in a nematic host into 2D square crystal-like 
structures was reported [29] . This has been achieved by 
mixing quadrupolar colloidal particles with distinct an- 
choring conditions, namely homeotropic and degenerate 
planar. The homeotropic anchoring particles nucleate 
Saturn-ring defects, and the planar anchoring particles 
antipodal boojums. In the linear regime each of these 
colloidal particles induce long-range distortions of the ne- 
matic director with quadrupolar symmetry (44|, and as 
a consequence the long-range pair interaction is of the 
quadrupolar type. In this paper we have studied numeri- 
cally the nematic-mediated interaction between two such 
particles. For simplicity we have assumed spherical par- 
ticles of equal radii. 

We have confirmed that the S and P particles have 
quadrupole moments wit h op posite signs, QsQp < 0. 
The asymptotic analysis [23| suggests that S particles 
have a positive quadrupole moment, and P particles a 
negative one. At large distances the interaction poten- 
tial F^{d,6) behaves as d~ a with a = 4.8 ± 0.2 which 
is in excellent agreement with the asymptotic result d~ 5 
for quadrupolar interactions. We have also confirmed nu- 
merically the experimental values for the product of the 
quadrupole moments QsQp = — 0.18i? 6 . Note that if the 
confining surfaces of the experimental setup (29j had an 
influence on the pair interaction, there would be a devia- 
tion from the power law behaviour that would eventually 
(for very strong confinement) turn into an exponential 
decay. Since in our calculations the finite size effects are 
negligible, our results indicate that in the experiment of 
Ognysta et al. [29| the confining surfaces played a negligi- 
ble role. For distances d < 3i?, the asymptotic quadrupo- 
lar result Eq.© is qualitatively wrong, see Fig. [5J At 



By contrast to the interaction of particles with iden- 
tical anchoring conditions (P-P, or S-S), the S-P inter- 
action potential F^ 2 \d,9) exhibits two minima (in the 
first quadrant), when the particles are oriented at 9 = 
or 9 = tt/2 relative to the far- field director, and sepa- 
rated by <f|| ~ 2.16R and d± ~ 2.25R, see Fig. [U These 
values are in excellent agreement with those reported ex- 
perimentally [H, d** p ~ 2.16i? and d e * p ~ 2.28R [Hj]. 
We have also found a strong repulsion at short distances 
and at any orientation 9. These two features of the pair 
interaction potential are pre-requisites for the assembly 
of colloidal lattices with rectangular unit cells. 

We have calculated the pair clastic force, see Fig. [3] 
The polar component Fg is always zero at 9 = 0, tt/2, 
and acquires rather large values at oblique orientations 9, 
provided that the inter-particle distance is not too large 
(d < 3i?, see Fig. [3] right panel) . The sign of Fg depends 
on 9. Fg < for 9 ^ 9 c (d), and it is positive other- 
wise. 9 c (d) is the angular position of the cusp in the 
F^ (9) profile, see Fig. SJ Therefore, there is a strong 
clastic force which aligns a pair of colloidal particles ei- 
ther parallel (9 ^ 9 c (d)), or perpendicular (9 ^ 9 c (d)) to 
the far-field director. Once the particles are aligned along 
one of these directions, they move under the action of the 
radial force F r towards the global d = dn,6 = 0, or the 
local d = d±,9 = tt/2 equilibrium positions. The pres- 
ence of a strong aligning force Fg at small-to- intermediate 
distances is in sharp contrast with the P-P interaction 
profile [22j , where the polar component Fg is always two 
orders of magnitudes smaller than the radial force F r . 

Most theoretical studies of the interaction between col- 
loidal particles in nematic liquid crystals focus on parti- 
cles with the same size. Although this may seem a rather 
strong simplification, it is justified by the current experi- 
mental ability of producing monodispersed colloidal par- 
ticles. This seems to be the best choice for the assembly 
of photonic crystals, but the effect of polydispersity on 
the stability of such crystals remains an important open 
question. 
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